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^_l ' Abstract 

We present a Friedmann-Robertson- Walker quantum cosmological model in the presence of Chaplygin 
gas and perfect fluid for early and late time epoches. In this work, we consider perfect fluid as an effective 
potential and apply Schutz's variational formalism to the Chaplygin gas which recovers the notion of time. 
These give rise to Schrodinger-Wheeler-DeWitt equation for the scale factor. We use the eigenfunctions 

j^i in order to construct wave packets and study the time dependent behavior of the expectation value of the 

scale factor using the many-worlds interpretation of quantum mechanics. We show that contrary to the 

^rv' classical case, the expectation value of the scale factor avoids singularity at quantum level. Moreover, this 

model predicts that the expansion of Universe is accelerating for the late times. 
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0\ ■ 1 Introduction 



Supernova la (SNIa) observations show that the expansion of the Universe is accelerating [1], contrary to 
Friedmann-Robertson- Walker (FRW) cosmological models, with non-relativistic matter and radiation. Also 



cosmic microwave background radiation (CMBR) data [2, 3] is suggesting that the expansion of our Universe 
seems to be in an accelerated state which is referred to "dark energy" effect [4]. Cosmological constant, A, as 
the vacuum energy can be responsible for this evolution by providing a negative pressure [5, 6]. Unfortunately, 
the observed value of A is 120 orders of magnitude smaller than the one computed from field theory models 
[5, 6]. Quintessence is an alternative to consider a dynamical vacuum energy [7], involving one or two scalar 
fields, some with potentials justified from supergravity theories [8]. However, the fine-tuning problem of these 
models which arises from cosmic coincidence issue has no satisfactory solution. 

The Chaplygin gas model is an interesting proposal [9], describing a transition from a Universe filled with 
dust-like matter to an accelerating expanding stage. This model was later generalized in Ref. [9, 10]. The 
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generalized Chaplygin gas model is described by a perfect fluid obeying an exotic equation of state [9, 10] 

where A is a positive constant and < a < 1. The standard Chaplygin gas [9] corresponds to a — 1. Some 
publications [11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29] and reviews [23, 24] which 
studied the Chaplygin gas cosmological models have already appeared in the literature. The Chaplygin gas 
can be obtained from the string Nambu-Goto action in the light cone coordinate [26]. Since the application 
of string theory in principle is in very high energy when the quantum effects is important in early Universe, a 
quantum cosmological study of the Chaplygin gas is also well founded. 

Recently, Quantum mechanical description of a FRW model with a generalized Chaplygin gas has been 
discussed in Rcf. [30] in order to retrieve explicit mathematical expressions for the different quantum mechanical 
states and determine the transition probabilities towards an accelerated stage. Moreover, quantization of FRW 
model in the presence of Chaplygin gas has been discussed in Ref. [31]. There, we have considered matter as 
the Chaplygin gas and discussed the early time behavior of expectation value of the scale factor through the 
application of Schutz's formalism. In this paper, aside from the Chaplygin gas which is coupled to gravity 
and has an advantage of furnishing a variable connected to matter which can be identified with time, we also 
include the perfect fluid in this scenario and investigate the analytical solutions in both early and late time 
Universes. Schutz's formalism [32, 33] gives dynamics to the matter degrees of freedom in interaction with 
the gravitational field. Using proper canonical transformations, at least one conjugate momentum operator 
associated with matter appears linearly in the action integral. Therefore, a Schrodinger-likc equation can be 
obtained with the matter variable playing the role of time. The application of Schutz's formalism in Stephani 
and FRW perfect fluid cosmological models has been discussed in Refs. [34, 35, 36]. Note that our approach in 
principle is different from Moncrat et. al [37], where they correspond the dynamical variable to the perfect fluid 
instead of Chaplygin gas and resort to the numerical methods to obtain the time evolution of an initial wave 
packet. Note that there is considerable evidence that the early Universe is dominated by radiation. Therefore, 
a natural setting for quantum cosmology is the one where radiation has the predominant role [38]. On the 
other hand the Chaplygin gas is dominated by the non-relativistic matter at early times (See the following 



section). This seems to be in contradiction with our knowledge of baby Universe. According to [39], inflation 
can be accommodated within the generalized Chaplygin gas scenario. Hence, the way adopted to avoid this 
inconsistency is that the radiation dominated phase is followed by Chaplygin dominated period so that we have 
the so-called Chaplygin inflation [39]. Also, it would be more suitable to consider the field theory representation 
of the Chaplygin gas [40] to describe quantum cosmology. In this way, the Chaplygin gas can be viewed as 
a modification of gravity as was first pointed out in [40]. Also, the authors of [41] have recently shown that 
the Chaplygin gas model has a geometrical explanation within the context of brane world theory for any a. 

Consequently in these models the equation 

i 
B 



A 



(2) 



a(i) 3 ( Q+1 ) 

is a consequence of stress-energy conservation for a scalar field on the brane [40] , or conservation of induced 
dark matter on the brane [41, 42]. Here, a(t) is the scale-factor of the universe and B is a positive integration 
constant. Therefore, it is evocative to view the contribution of Chaplygin gas to the stress-energy tensor as a 
brane induced modification of gravity. In this article, we used the fluid description for the Chaplygin gas and 
for Lagrangian formalism, the corresponding pressure. Consequently, if we rely on the model described in [41], 
we will have covariance in our model. 

The paper is organized as follows. In Sec. 2, the quantum cosmological model with a Chaplygin gas, as a 
portion of the matter content is constructed in Schutz's formalism [32, 33] for early and late time Universes. 
Then the Schrodinger-Wheeler-DeWitt (SWD) equation in minisuperspace is obtained to quantize the model 
under the action of a perfect fluid effective potential. The wave function depends on the scale factor a and on 
the canonical variable associated to the Chaplygin gas, which in the Schutz's variational formalism plays the 
role of time T. We separate the wave function into two parts, one depends only on the scale factor and the 
other depends on the time. The time dependent part of the solution is e lEt , where E is the energy. In Sec. 3, 
we construct wave packets from the eigenfunctions and compute the time-dependent expectation values of the 
scale factor to investigate the existence of singularities at quantum level. Moreover, we present some analytical 
solutions in both early and late time epoches. In Sec. 4, we present our conclusions. 



2 The Model 

The action for gravity plus Chaplygin gas in Schutz's formalism is written as 

S= f d 4 xy/=gR + 2 f d i x\fhh ab K ab + f d i x^/^p f + f d 4 xy/=g p c , (3) 

JM JdM JM JM 

here, K ab is the extrinsic curvature and h a \, is the induced metric over the three-dimensional spatial hypersur- 
face, which is the boundary dM of the four dimensional manifold M. We choose units such that the factor 
167rG becomes equal to one. pf and p c denote Chaplygin gas and perfect fluid pressure, respectively. Note 
that, according to [32] the above action is equivalent to the usual Hawking-Ellis formalism for perfect fluid 
description [43] . Perfect fluid satisfies the barotropic equation of state 

Pf = wpf, w < 1. (4) 

The first two terms were first obtained in [44] and the last term of (3) represents the matter contribution to the 

total action. In Schutz's formalism [32, 33] the fluid's four- velocity can be expressed in terms of five potentials 

$, C, 0, and S 

u v = -{$, v + 0, v + 6S, v ) (5) 

A* 

where \x is the specific enthalpy. S is the specific entropy and the potentials £ and (3 are connected with 
rotation which are absent of models in the Friedmann-Robertson- Walker (FRW) type. The variables <& and 9 
have no clear physical meaning. The four-velocity also satisfies the normalization condition 

u v u u = — 1. (6) 

The FRW metric 

ds 2 = -N 2 (t)dt 2 + a 2 {t)g t3 dx l dx\ (7) 

can be inserted in the action (3), where N(t) is the lapse function and gij is the metric on the constant-curvature 
spatial section. Following the thermodynamic description of Ref. [45], the basic thermodynamic relations take 
the form 

Pc = p [l + n], h = l + U + p c /p , (8) 



TdS 



dn + Pc d(i/po), 


l + a 


(i + n) 1+Q + -^ 

L Po 



(9) 



It then follows that to within a factor 



(i + n)-° 

l + a ' 

s = (i + n) 1+a - 



A 



(10) 
(11) 



Therefore, the equation of state takes the form 

Pc = -A 



A I * S-V" 



(12) 



The particle number density and energy density are, respectively, 



Pc = 



Po 






-1 

l + a 



(13) 
(14) 



where /i = ($ + 8S)/N. After dropping the surface terms, the final reduced action takes the form 



S = J dt\ - 6 + QkNa - Na A p f - Na 3 A 



1 I (<fr + 0S) — 



(15) 



The reduced action may be further simplified using canonical methods [45] , resulting in the supcr-Hamiltonian 



6ka + a 3 Pf + (Spl +a + Aa 3< - 1+a ^ 



n = ~24a~ {>1 ""^" r ' f l " i ' ,J ' 



OC 



(16) 



where p a = —12da/N and p$ = ^z- . However, an analytical quantum mechanical treatment of this FRW 
minisupcrspace with the above Hamiltonian does not seem feasible. Therefore, it requires some approximation. 
We study the Chaplygin gas expression in early and late times limits, namely for small scale factors Sp$ a 3> 
^ fl 3(i+a) [30^ 31] and large scale factors Sp]^ a <C Aa 3i ^ 1+a \ separately. So for early Universe, we can use the 
following expansion 



(Sp 



i+« + Aa m+*)) — 



S!+«p$ 



1 Aa 3(a+1 *> 1 1 



1 + a Spl +a 21 + a \l + a 

Hence, up to the leading order, the supcr-Hamiltonian takes the form 

H = -— - 6ka + a 3 p f + S^p®. 
24a J 



A 2 



s 2 P t 



2(l+a) 



,6(a+l) 



.(17) 



(18) 



The following additional canonical transformations 



T=-(l + a)p^ 1 S l +"p S) p T = Si+op®, 



(19) 



and use of the explicit form of the energy density of the perfect fluid pf = 3 , 1+ > , simplify the super- 



Hamiltonian to 



V 2 
H = - — -6ka + Ba~ 
24a 



■Pt, 



(20) 



where B is a constant and the momentum px is the only remaining canonical variable associated with matter. 
It appears linearly in the super-Hamiltonian. The parameter k defines the curvature of the spatial section, 
taking the values 0, 1, —1 for a flat, positive-curvature or negative-curvature Universe, respectively. 

The classical dynamics is governed by the Hamilton equations, derived from Eq. (20) and Poisson brackets 



MS 



' a= {a,7VH} = -^, 

Pa= {pa,NH} = -^p 2 a + 6Nk + 'iwNBa 

t = {T, NH} = N , 

p T = {p T ,NH} = 0. 



(21) 



We also have the constraint equation Ti = 0. Choosing the gauge N = 1, we have the following solutions for 
the system 



T = t, 



Pt = const., 

a = -7T- -7T- -rwBa- 3w - 2 , 
2a 2a 4 

= -6aa 2 - 6ka + Ba^ w + p T . 



(22) 
(23) 
(24) 
(25) 



The WD equation in minisuperspace can be obtained by imposing the standard quantization conditions on the 

Br) 
canonical momenta (p a = —i-M-, Pt = ~i~pfp ) an d demanding that the super-Hamiltonian operator annihilate 



the wave function (h = 1) 



d 2 ^ d'fy 

—-^- - (lUka 2 - 24Ba 1 - 3w )y - i24a— = 0. 
oa z at 



(26) 



In this equation according to (22), T — t corresponds to the time coordinate. As discussed in [46, 47], in order 
for the Hamiltonian operator H to be self-adjoint the inner product of any two wave functions $ and ^ must 
take the form 

/■OO 

($,*) = / a$*t>da 7 (27) 

Jo 

On the other hand, the wave functions should satisfy the following boundary conditions 

&f>(a,t) 



*(0,t) = or 



da 



0. 



a=0 



The SWD equation (26) can be solved by separation of variables as follows 



(28) 



#z,t) = e l ^{a), 



(29) 



where the a dependent part of the wave function tp(a) satisfies 



-tp"(a) + (144/ca 2 - 24Ba 1 - 3w )iP(a) = 24Eaijj(a), 



(30) 



and the prime means derivative with respect to a. 

Now, we consider late time Universe when Sp$~ a <C ^4a 3 ( 1+a ) . Using the expression 

1 Spi +a , 1 1 



(Spl +a + Aa 3< - 1+a ^ « AT^a 3 



l + aAa^ a+1 ) 21 + all + a J A 2 a^ a+1 ) 



1 A gpg +a) 



up to the first order, the supcr-Hamiltonian (16) takes the form 



A 



3 i A „-3an 1+a 



H = -— - Qka + aTpf + A— a 6 + a~ ia Sp. 

24a ' 1 + a 



$ ! 



(31) 



(32) 



The following additional canonical transformations 



(!+«)„ 



A— 



T = -(1 + a)A-— P$ ^> s , PT = —— Sp. 

1 + a 



3 $ 7 



(33) 



simplify the super-Hamiltonian to 



n 2 i 

W = --£2- - 6/ca + Ba" 3tu + A~ a 3 + a~ 3 > T . 
24a 



(34) 



The classical dynamics is governed by the Hamilton equations, derived from Eq. (20) and Poisson brackets as 



(35) 



a = { a ,NH} = -^., 

Pa= { Pa ,NH} = - J ^p z a + 6Nk + 3wNBa- 

SNA^a 2 + 3aN a" 3 "" 1 ^, 
t= {T,NH} = Na- 3a , 
p T = { PT ,NH} = 0. 

We also have the constraint equation H = 0. Choosing the gauge N — a 3a , we have the following solutions for 
the system 



t, 



PT — const. 



3a _ « 

la I 



- -wBa 6a - 3w - 2 + -ATha 6a+1 - -ap T a 3a - 2 , 



= 



-6a- 6a+1 a 2 - (Ska + B a - 3w + A^ a 3 + a~ 3a p T . 



(36) 

(37) 
(38) 
(39) 



It is important to note that these equations predict an accelerating Universe for late times. For large values of 
the scale factor we can simplify the above equations and find the acceleration parameter 



q = — = da — 1, 



(40) 



which is positive for a > 1/3. Now, imposing the standard quantization conditions on the canonical momenta 
and demanding that the supcr-Hamiltonian operator annihilates the wave function, we are led to SWD equation 
in minisupcrspacc (K = 1) 



a 2 * 



- (144fca 2 - 2<±Ba x - 3w - 24.4.— a 4 )* _ i24a 



1-3^* 



= 0. 



(41) 



da 2 K dt 

Here, according to (36), T — t corresponds to the time coordinate. Demanding that the Hamiltonian operator 
H to be self-adjoint, the inner product of any two wave functions $ and $ must take the form [46, 47] 



(*,*) 



a l-3a $ *^ rfa _ 



(42) 



The SWD equation (41) can be solved by separation of variables as follows 



^(o,t) = e Ul V(o), 



(43) 



where the a dependent part of the wave function ip{a) satisfies 



-ip"(a) + (lUka 2 - 2ABa 1 - Zw - 24A^ a 4 \ ip(a) = 2AEa 1 - 3a ip(a), 



(44) 



and the prime means derivative with respect to a. Note that effective Chaplygin gas term (24A 1 +° ) plays the 
role of a positive cosmological constant. In particular, when a — 1/3 and w — 1/3, this equation reduces to 
the FRW model with positive cosmological constant and radiation which has been studied in Ref. [48] . 

3 Results 

In this Section we first study the issue of singularity avoidance in quantum cosmology in the early Universe 
and then present some analytical solutions in both early and late Universes. 

For k = the time- independent Wheeler-DeWitt equation (30), in the dust dominated Universe (w = 0), 
reduces to 

ip" + 24(E + B)aiP = 0. (45) 

The above equation has the following general time-dependent solutions under the form of Bessel functions 



V E = e* Et ^ 



C\ J 



(V96E 7 A TA /V96S 7 
i — ^ — a 2 +c 2 Yi — - — a 



3 V 3 



3 V 3 



(46) 



where E' = E + B. Now, the wave packets can be constructed by superimposing these solutions to obtain 
physically allowed wave functions. The general structure of these wave packets are 



/■OO 

*(a,t)=/ A(E')^' E (a,t)dE'. 

Jo 



(47) 



We choose C2 = for satisfying the first boundary condition (28). Defining r — ^ , simple analytical 
expressions for the wave packet can be found by choosing A(E') to be a quasi-gaussian function 



/■OO 

Jo 



(48) 



where v = -| and 7 is an arbitrary positive constant. The above integral is known [49], and the wave packet 



takes the form 



*(o,t) = a- 



(-2Z): 



(49) 



Figure 1: The time behavior of the expected value for the scale factor (a)(i) (solid line) and the classical scale 
factor a(t) (dashed line) for dust dominated Universe (w = 0) and flat space time (k = 0). 

where Z = 7— i -^t. Now, we can verify what these quantum models predict for the behavior of the scale factor 
of the Universe. By adopting the many- worlds interpretation [50, 51], and with regards to the inner product 
relation (27), the expectation value of the scale factor 



< a> (t) = 



J °° a^{a,t)*^{a,t)da 



is easily computed, leading to 



< a > (t) ex 



:t z + 1 



(50) 



(51) 



,(32) 2 7 2 
These solutions represent a no singular Universe which goes asymptotically over to the corresponding flat 

classical model for dust (w = 0) dominated epoch (22-25) (Fig. 1) 



a(t)o,t 2 '\ 



(52) 



In the case k = 1 and w = the time-independent Wheeler-DeWitt equation (30) reduces to 

-tp"(a) + {-2AE'a + 144a 2 ) r/>(a) = 0. 
Defining new variable x — 12a — E' we find 



d 2 4> 
dx 2 



E' 2 x 2 



ip{a) = 0. 



(53) 



(54) 



144 144_ 

Equation (54) is similar to the time-independent Schrodingcr equation for a simple harmonic oscillator with 
unit mass and energy A 

(55) 



d 2 ip 



-^ + [-2A + wV]V(aO = 0, 



10 



where 2A = _E' 2 /144 and w — 1/12. Therefore, the allowed values of A are w(n + 1/2) and the possible values 
of E' are 



E' n = y/\2(2n+l) , n = 0,l,2,... . (56) 

therefore, the stationary solutions are 

* n (o, i) = e i£ "Vn (12a - E' n ) , (57) 

where 



x 



V n (x) = H n {^=)e-*^, (58) 



and i?„ arc Hcrmite polynomials. The wave functions (57) are similar to the stationary quantum wormholcs as 
defined in [52]. However, neither of the boundary conditions (28) can be satisfied by the these wave functions. 
In k = — 1 and w = case, equation (30) reduces to 

i>"(a) + {2AE'a + 144a 2 ) r/>(a) = 0, (59) 

where the solutions are 

I 48 '4 y 12 y 48 U y 12 / J 

where M K} \ and W k ,a are Whittaker functions. The Whittaker functions do not automatically vanish at a = 0. 
Therefore, we need to take both C\ ^ and C2 7^ to satisfy ^(0, t) = 0. 
For w = —1/3, the SWD equation (30) can be written as 

-4>"(a) + 24(6fc - B)a 2 iP(a) = 24Eaip(a), (61) 

which as before, has the solutions in the form of Simple Harmonic Oscillator (57) with discrete spectrum or 
Whittaker function (60) for positive or negative value of (6/c — B), respectively. 
For k = and w = 1/3 (radiation), the WD equation (30) reduces to 

-ip"{a) - 24BV(a) = 24EaiP(a), (62) 

which can be rewritten as 

V>"(a) + 24E (a + ^ J ip(a) = 0, (63) 

11 



D 



by taking x = a + % we have 



dx 2 



ip{x) + 24ExiP(x) = 0, 



(64) 



which is the Airy's differential equation. We solve this equation for E > and E < 0, separately. 

For E > 0, this equation has two solutions as Ai [— (24E) 1 ^ 3 x\ and Bi [— (24E)) 1 ^ 3 x\ . First one is ex- 
ponentially decreasing function of x and the second one grows exponentially and is physically unacceptable. 
Therefore, the solution is 



tp(a) = Ai 



-(2425)) 1 " [a+- 



D 



E 



We choose the first boundary condition (28), which leads to 



Ai 



-(24E) 



1/3 



£ 

£ 



0. 



(65) 



(66) 



Airy's function Ai(x) has infinitely many negative zeros z n — — a„, where a n > 0, therefore, the energy levels 

quantize and take the values 

{24 1 / 3 bV /2 
E n = *=-?- • (67) 

y a n J 

The time-dependent eigenfunctions take the form 

B 



*„(a,t) = e^ t Ai 



-{2AE n 



.1/3 



E,. 



(68) 



For E < 0, this equation has also two solutions as Ai [(24|i?|) 1 / 3 a;] and Bi [(24|S|)) 1 / 3 a:] . Since the second 
one grows exponentially and is physically unacceptable, the solution is 

B 



(24|£|)) 1/3 U 



tp(a) = Ai 



We choose the first boundary condition (28), which leads to 



Ai 



vl/3 



B 



-(2A\E\y'* — 



\E\ 



= 0, 



(69) 



(70) 



therefore, the energy levels quantize and take the values 



E„ 



The time-dependent eigenfunctions take the form 



2A l ' 3 B 



.3/2 



(71) 



* n {a,t) = e l ^ t hi 



(24KI) 1 / 3 [a 



B 



(72) 
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Figure 2: Plot of the wave function (ip(a)) for B — 1 and n — 8, showing the oscillatory behavior for the small 
values of the scale factor and exponential damping for the large values of the scale factor. 

It is important to note that Airy's function Ai(x) has an oscillatory behavior for x < (a < rjgr-r ) whiles 

\E n \ 

for x > (a > jttt) decreases monotonically and is an exponentially damped function for large x (Fig. 2). 

\E n \ 

Therefore, the solutions (72) show a classical behavior for small a and a quantum behavior for large a. This is 
contrary to usually expected results for previous case. In fact detecting quantum gravitational effects in large 
Universes is noticeable which has been also observed in FRW, Stephani, and Kaluza-Klein models [53, 35, 54]. 



In k = 1 and w = 1/3 (radiation) case, the WD equation (30) reduces to 



-i>"(a) + (144a 2 - 2AB)^{a) = 24£aV>(a). 



(73) 



The above equation can be written as 

-V>"(a) + 144 
by taking x = a — jj we have 



-!Hi) 2 



V = o, 



(74) 



—r-zi>(x) + UAx 2 xP(x) = (E 2 + 24BU(x). 
dx z 



(75) 



This equation is identical to the time-independent Schrodinger equation for a simple harmonic oscillator with 
unit mass and energy A 



d 2 ip(x) 



+ uj 2 x 2 tp(x) = 2\tp(x), 



(76) 



where 2A = (E 2 + 24B) and lo 2 = 144. Therefore, the allowed values of A are u)(n + 1/2) and the possible 
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values of E are 



E n = v/6(n + 1/2) - 24B , n = 0,l,2,... , (77) 

therefore, the stationary solutions are 

^(a,t) = e i£ »Vn(a-|), (78) 



^„(x)=i7„^12)2xj e - 3a; , (79) 

where H n are Hermite polynomials. However, neither of the boundary conditions (28) can be satisfied by these 
wave functions. 

Now, we present some analytical solutions for the late time Universe. For flat space time (k — 0), dust 
epoch (w = 0), and standard Chaplygin gas (a = 1), equation (44) reduces to 

-ip"(a) + (-24Ba - 2AA^a A ) ip(a) = 24£a~ 2 V(a), (80) 

where the solutions are 



*(0,t) - e *^ e «vlA»(i+«.)a« l-iVI=58B 



X 



Ci [/ [ - (-2iV6BA W^> - Vl - 96£ + 3V 1 - -\/I - 96£, -4iW -A^r+^T a 3 



--|V^96£ / A .J2 A _^_ ■, 



+ c 2 L-av— _ x _4i A /rA5P+=y d . (81) 

i(2i\/6BA 2 ( 1 +°) +V1-96B-3) \ » O y J 

Here £/(a, 6, c) is the confluent hypergeometric function and L^{x) is the generalized Laguerre polynomial. We 
need to take both C\ ^ and C 2 ^ to satisfy \I>(0, t) = 0. 

For flat space time (k = 0), stiff matter (w = 1), and standard Chaplygin gas (a — 1), equation (44) reduces 
to 

-ip"(a) + (-24Ba- 2 - 24A^ a 4 ) tp(a) = 24£a" 2 V(a), (82) 

with the solutions as 

*(«•*) - etEt { C ^ aJ -i^M^^( 2 \[l AJuUa3 ) 



+ C ' ^ J W -96 ( B +E) + l [ 2 \l3 A ~ a3 ) )■ ( 83 ) 

Here again, we have C\ 7^ and Ci =/= in order to satisfy the first boundary condition (28). 
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4 Conclusions 

In this work we have investigated minisuperspace FRW quantum cosmological models with Chaplygin gas and 
perfect fluid as the matter content in early and late times. The use of Schutz's formalism for the Chaplygin 
gas allowed us to obtain SWD equations with the perfect fluid's effective potential. We have obtained eigen- 
functions and therefore acceptable wave packets were constructed by appropriate linear combination of these 
cigenfunctions. The time evolution of the expectation value of the scale factor has been determined in the 
spirit of the many- worlds interpretation of quantum cosmology. We have showed that contrary to the classical 
case, the expectation value of the scale factor avoids singularity at the quantum level. Moreover, this model 
predicts an accelerated Universe for late times. 
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